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Aims. The main aim is to derive reliable mass-loss rates and circumstellar SiO abundances for a sample of 40 S-type AGB stars based 
on new multi-transitional CO and SiO radio line observations. In addition, the results are compared to previous results for M-type 
AGB stars and carbon stars to look for trends with chemical type. 

Methods. The circumstellar envelopes are assumed to be spherically symmetric and formed by a constant mass-loss rate. The mass- 
loss rates are estimated from fitting the CO observations using a non-local, non-UTE radiative transfer code based on the Monte-Carlo 
method. In the excitation analysis, the energy balance equation is solved self-consistently simultaneously as the radiative transfer and 
the temperature structure of the gas is derived. Effects of dust grains are also included in the molecular excitation analysis. Once the 
physical properties of the circumstellar envelopes are determined, the same radiative transfer code is used to model the observed SiO 
lines in order to derive circumstellar abundances and the sizes of the SiO line-emitting regions. 

Results. We have estimated mass-loss rates of 40 S-type AGB stars and find that the derived mass-loss rates have a distribution that 
resembles those previously derived for similar samples of M-type AGB stars and carbon stars. The estimated mass-loss rates also 
correlate well with the corresponding expansion velocity of the envelope, in accordance with results for M-type AGB stars and carbon 
stars. In all, this indicates that the mass loss is driven by the same mechanism in all three chemical types of AGB stars. In addition, 
we have estimated the circumstellar fractional abundance of SiO relative to H2 in 26 of the sample S-type AGB stars. The derived SiO 
abundances are, on average, about an order of magnitude higher than predicted by stellar atmosphere thermal equilibrium chemistry, 
indicating that non-equilibrium chemical processes determines the abundance of SiO in the circumstellar envelope. Moreover, a 
comparison with the results for M-type AGB stars and carbon stars show that for a certain mass-loss rate, the circumstellar SiO 
abundance seems independent (although with a large scatter) of the C/O-ratio. 

Conclusions. In our comparison of S-type AGB stars with carbon stars and M-type AGB stars, we find no large differences in 
circumstellar physical properties or SiO abundances depending on the chemical type of the star 

Key words. Stars: AGB and post-AGB - Stars: abundances - Stars: carbon - Stars: late-type - Stars: mass-loss 



1. Introduction 

The final evolutionary stage of low- to intermediate-mass stars, 
as they ascend the asymptotic giant branch (AGB), is character- 
ized by an intense mass loss. The stellar wind builds up grad- 
ually and creates a circumstellar envelope (CSE), carrying gas 
and dust from the star into the interstellar medium. The molec- 
ular setup and grain types in CSEs are to a large extent deter- 
mined by the C/O-ratio in the photosphere of the central star 
The AGB stars are normally divided into two distinct spectral 
types: M-type stars, with C/0<1 and carbon stars, with C/0>1. 
Abundance analysis has also revealed stars with photospheric 
C/O-ratios close to unity (w ithin approximately 5 %), so called 
S-type AGB stars ( IScalo & R oss 197Q). The spectra of S-type 
AGB stars are dominated by ZrO bands (as opposed to spectra of 
M-type AGB stars which are dominated by TiO bands) indicat- 
ing that the S-type AGB stars have an enhancement in elements 
formed through the slow neutron capture process. 

Due to having a C/O-ratio close to unity, it is tempt- 
ing to identify S-type AGB stars with a brief transitional 
phase as the star evolves from an oxygen-rich M-type AGB 
star into a carbon star. Dredge-up of carbon from He-shell 



burning would change the spectral type of the star sequen- 
tially: M-MS-S-SC-C. As possible transition objects, S-type 
AGB stars might very well help to achieve a deeper un- 
derstanding of the chemical evolution as a star ascends the 
AGB, as well as shed light on the mass-loss mechanism(s), 
which is (are) not yet fully understood in detail. Several sur- 
veys of CO emission from their CSEs have b e en performed 
(Biedng & L atteii Il994t ISahai & Liechtil 119951; rBieging et al.l 
1998,; ,Groenewegen & de Joiii 119981; iJorissen & Knapplll998l; 



Ramstedt et al. I2OO6I) . Circumstellar molecular line 
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from molecules other than CO has previously been searched for, 
and detected, only in a handful of objects and onl y for HCN and 
SiO dBieging & Latteilll994l;lBieging et al.lll998l) . 



The aim of this work is to thoroughly investigate the circum- 
stellar physical and chemical properties of a sample of S-type 
AGB stars. The physical properties of the CSEs, in particular 
the mass-loss rates that produced them, are determined from the 
CO data using detailed, non-LTE, radiative transfer modelling 
which self-consistently calculates also the gas kinetic tempera- 
ture. Detaile d studies of the mass-loss properties of sampl es of 
carbon stars dSchoier & Olofssonll200ll ; Schoier et al1l2002l). M- 
type AGB stars (IOlofssonetalJl2002t iGonzalez Delgado et al.l 
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l2003h . and S-type AGB stars dRamstedt et al.ll2006l) have previ- 
ously been performed. 

Lately, the mechanisms driving mass loss in M- 
type AGB stars have been much deb ated (IWoitkd 120061: 
lH5fner & AndersenI 120071: iHofneij |2008|) . and for the S-type 
AGB stars this is a matter of a long-standin g debate (e.g., 
IWillems & de Jong|ll988l:ISahai & Liechtilll995h . The supposed 
lack of any free oxygen and/or free carbon to form dust (due to 
a C/O-ratio close to unity) has lead to the suggestion that the 
S-type AGB stars would not be able to drive a wind as efficiently 
as the M -type AGB star s and t he carbon stars. However, the 
results of iR amstedt et al.l ^2006) show no apparent differences 
between the mass-loss rate distribution of the S-type AGB stars 
and those found for carbon and M-type AGB stars. 

Once the physical properties are known, from the CO 
analysis, they can be used to estimate abundances of other 
molecul es in the CSE . The analysis of the S-type AGB 
stars in iRamstedt et alj (|2006|) was based on observations of 
CO(7 = 3— >2) data gathered at the APEX telescope supple- 
mented with data collec ted from the lit erature. When examin- 
ing the published data, IRamstedt et alj (|20()6) found a rather 
large scatter in the reported line intensities [especially for 
CO(7 = 2— > 1) data] for individual objects, even when observed 
with same telescope. Since then, we have obtai ned more obser- 
vation al data, which warrants a re-analysis of the lRamstedt et alj 
(120061) work, and the results are reported here. 

We have, as a first step after the CO observations, searched 
for circumstellar SiO radio line emission in several rotational 
transitions. Using the physical structure of the CSEs derived 
from the CO modelling, we estimate the abundance of SiO in 
the same sample of S-type AGB stars based on a detailed ex- 
citation analysis. There exists strong evidence that the circum- 
stellar SiO emission carries information on the region where 
the m ass loss is initiated and where the dust formation takes 
plac e (Reid &Moranll 19811: iReid & Menteiill 19971: ISchoier et al] 
2006a) making it a p articularly interesting molecule to study. 
Schoier et alj ( l2006bl) modelled circumstellar SiO line observa- 



Table 1. The sample with variable types, periods P, luminosities 
L*, and distances D. 



tions from a sample of carbon s tars and when comparing this to 
a similar analysis performed bv lGonzalez Delgado et al.l ( l2003l) 
for a large sample of M-type AGB stars, they found no ap- 
parent distinction between their circumstellar SiO abundance 
distributions. For the carbon stars, the derived abundances are 
several orders of magnitude higher than expected from thermal 
equilibrium stellar atmosphere chemistry. A possible explana- 
tion for the high SiO abundances derived for the carbon stars 
is the influence of a sh ock chemistry in the inner part of the 
wind dCherchneffl |2006|) . With the analysis performed in this 
work, the S-type AGB stars can be added to this comparison 
(see Sect.|63]). 

The sample of S-type AGB stars and the observational data 
are presented in Sects |2| and |3] The radiative transfer modelling 
is described in Sect. H] The results are given and discussed in 
Sects |5] and |6] respectively. Finally, the conclusions are given in 
Sect. Ill 



2. The sample 

iJorissen & Kn arip (19 9j|) presen ts a sample of 124 S-type stars 
from the list o f Chen et al.l (Il995h providing cross identifications 
between the General Catalogue of Galactic S stars, the IRAS 
Point Source Catalogue (PSC), and the Guide Star Catalogue. 
Only stars having flux densities of good quality in the 12, 
25, and 60 fim bands in the IRAS PSC were retained in their 
sample. We have chosen the stars in the sample of Jorissen 
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Lb 




4000 


800 


TCam 


M 


374 


5600 


540 


S Cas 


M 


611 


8000 


440 


V365 Cas 


SRb 


136 


4000 


625 


WY Cas 


M 


477 


6700 


600 


AMCen 


Lb 




4000 


750 


TTCen 


M 


462 


6500 


880 


UYCen 


SR 


115 


4000 


590 


V386 Cep 
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& Knapp that have previously been detected in circumstellar 
CO emission. Five stars have b een added from the sample s of 



Sahai&Liechtil(ll995h[se- 



Stephensonl (1 1 984 1 1 99 



iGroenewegen &de.Tongl(II998h and 

lected from the S-type catalogs of .^.^^^^^„.,v,i„ ,> ,. 
and on e star was added from th e sample of iBieging & Latteil 
(11994 [selected from the Hst of I Jural (Il988h l. These six stars 
have also been detected previously in circumstellar CO emis- 
sion. Only S-type stars with Tc lines detected in their spectra and 
with detectable infrared excess are safely identified as intrinsic 
thermally pulsing AGB stars. S-type stars with no Tc lines and 
no infrared excess are most likely extrinsic S-stars, i.e., they are 
part of binary systems and their chemical peculiarities are due 
to mass transfer across the system. Jorissen & Knapp investi- 
gate the extrinsic/intrinsic characteristics of their sample and we 
have only retained stars which they classify as intrinsic. The stars 
from the Groenewegen & de Jong and Sahai & Liechti samples 
have all been detected in Tc and/or show infrared excess. We 
therefore believe that our sample only consists of intrinsic S- 



S. Ramstedt et al.: Circumstellar molecular line emission from S-type AGB stars: Mass-loss rates and SiO abundances 



3 



type AGB stars. Our sample of 40 stars is listed in Table [T] Our 
selection criteria might introduce a bias in the sample toward 
S-type stars that have higher mass-loss rates, and S-type stars 
without or with very little mass loss will be missed. The reason 
for this selection is to be able to detect line emission from other 
molecules than CO. 




300 600 900 1200 

Distance [pc] 



Fig. 1. The distance distribution of our sample of S-type AGB 
stars ( squares) overlayed with the distribution found in lJural 
(Il990h assuming a smooth surface distribution (solid line, see 
text in Sect. |2]) of S-type AGB stars in the solar neighbourhood. 
The dashed line shows a fit to the observational results. The er- 
rors are based on Poisson statistics. 



For the Miras (M), bolometric lu minosities are es t imate d 
using the period-luminosity relation of IWhitelock et al.1 (Il994h . 
and periods are taken from the General Catalogue of Variable 
Stars (Kholopov et al. 1999). For the semiregular (SR, SRa, and 
SRb) and iiTegular variables (Lb), a luminosity of 4000 was 
assumed in accordance with lOlofsson et all (|2002|) . Corrections 
for the interstellar extinct ion, based on their galac tic latitude, 
were applied for each star (iGroenewegen et anil992l) . Distances 
are derived by fitting the spectral energy distribution (SED) cal- 
culated with DUSTY (Sect.g^l), to observed fluxes (Sect.|T4li. 
The derived distances, and adopted periods and luminosities are 
presented in Table [T] together with the variability type. The de- 
rived distances from the SED modelling agree well with dis- 
tances from Hipparcos parallaxes when available (for 7 stars). 

Fig. [T] shows the distance distribution of our sample in 150 
pc bins. The distance distribution is overlayed with a distribution 
assuming that the stars are evenly distributed, that the surface 
density of carbon stars is 40 kpc"^ with a scale height of 200 pc, 
and that there are 1/3 a s many S-type stars as carbon s tars in 
the solar neighbourhood dWing & Yorkall977l:|juralll990h . A di- 
rect fit to our sample within 600 pc, gives a surface density of 
20kpc"- with a scale height of approximately 150pc (dashed 
line. Fig. [TJ. Considering the uncertainties involved in the dis- 
tance estimate, classification, etc., we find the results to be con- 
sistent and believe that our sample is representative of mass- 
losing S-type AGB stars and complete to about 600 pc. The 
sample of carbon stars studied by Schoier & Olofsson (2001) is 
representative of C-type AGB stars and complete to about 500 
pc, while the completeness of the M-type AGB star samples of 



Table 2. Energies of the upper levels, E^p, main-beam efficien- 
cies, ?7mb, and main-beam FWHM:s, 0,nb, at the frequencies of 
the observed CO and SiO lines. 
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0.74 


22 
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1) 


230.538 


16.6 


0.52 


11 


IRAM 
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345.796 


33.2 


0.56 
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JCMT 


CO{J = 


3 




2) 


345.796 


33.2 


0.70 


18 


APEX 


Si0(y = 
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— > 


1) 


86.847 


6.3 


0.65 


44 


OSO 


siO(y = 
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— > 


1) 


86.847 


6.3 


0.78 


29 


IRAM 


SiO(7 = 


5 


— » 


4) 


217.105 


31.3 


0.57 


12 


IRAM 


SiO(7 = 


6 


— » 


5) 


260.518 


43.8 


0.69 


21 


JCMT 


SiO(7 = 


8 


— > 


7) 


347.331 


75.0 


0.63 


14 


JCMT 


SiO(7 = 
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—> 


7) 


347.331 


75.0 


0.70 


18 


APEX 



lOlofsson et aP (l2002h and iGonzalez Delgado et alJ (l2003l) have 
not been thoroughly investigated. 

3. Observational data 

When modelling circumstellar molecular line emission it is es- 
sential to observe as many different transitions as possible of the 
molecule under study. The derived mass-loss rate (when assum- 
ing a spherically symmetric, smooth wind) should be considered 
as the average mass-loss rate that created the CSE probed by 
the observed line emission. Different transitions probe slightly 
different regions in the CSE depending on their excitation re- 
quirement and thus the average is taken over a larger part of the 
CSE if more transitions are obtained. The same is true when esti- 
mating the radial abundance distribution and, in particular, when 
estimating the size of the emitting region of a specific molecule. 
In this work, we have tried to obtain observations of at least three 
transitions of CO and SiO for as many of the sample sources as 
possible (Tables lAlM^ . 

3.1. Observations of circumstellar CO and SiO radio line 
emission 

The first analysis of CO radio li n e emi ssion from this sample 
was published in iRamstedt et al.l (|2006|) . There, data from the 
literature (mainly J-l—>0 and 2^1) was analysed together 
with / = 1 ^ data collected with the Onsala Space Observatory 
(OSO) 20 m telescopf] (January 2006), and 7 = 3 ^ 2 data of 1 8 
sources collected with the APEX 12 m telescopeQ (from August 
to October 2005). 

The CO data base has since then been substantially extended. 
New observations of / = 1 ^ and 2 — > 1 line emission were 
performed at the IRAM 30 m telescop^ using the AB SIS re- 

' The Onsala 20 m telescope is operated by the Swedish National 
Facility for Radio Astronomy, Onsala Space Observatory at Chalmers 
University of Technology, with support from the Swedish Research 
Council. 

- The Atacama Pathfinder Experiment (APEX) is a collaboration 
between the Max-Plank-Institute fiir Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. 

3 IRAM is supported by CNRS/INSU (France), the MFC (Germany), 
and the IGN (Spain). 
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ceiver combination observing both polarizations simultaneously 
in August 2006. Further observations of 7 = 3 — > 2 line emission 
were performed at APEX using APEX-2A duringthe autumn 
of 2006 and at the James Clerk Maxwell Telescop^ (JCMT) in 
April 2007. At the JCMT we used HARP due to problems with 
the B receiver As spectrometers, the VESPA autocorrelator at 
IRAM, the EFTS at APEX, and the ACSIS digital autocorrelator 
at the JCMT were used. 

New SiO observations of 7 = 2 — > 1 line emission were per- 
formed at OSO using an SIS receiver and the autocorrela- 
tor spectrometer in the low-resolution mode in January 2006. 
Observations of 7 = 2 — > 1 and 5 — > 4 line emission were per- 
formed at the IRAM 30 m telescope in August 2006, and fur- 
ther observations of 7 = 6 — > 5 and 8-^7 line emission were per- 
formed at JCMT during the autumn of 2006 and spring 2007. 
The observations at IRAM and the JCMT were obtained si- 
multaneously as the CO observations using the same receiver 
and spectrometer setups. During 2006 and 2007 observations of 
7 = 8 — » 7 line emission was performed at APEX using the same 
setup as for the CO observations (all these lines are in the v = 
state). Simultaneously as the 7 = 2^1, v = 0, observations at 
OSO, the 7 = 2— >1, v - 1 line (which is of maser origin) was 
observed. We observed SiO(7 = 2— »1, v=l) line emission in 24 
of our sample stars. For four of the stars, WY Cas, V386 Cep, 
TV Dra and R Gem, this is the first detection of SiO maser emis- 
sion. 

At OSO the observations were performed using the dual 
beam-switch mode. The observations at IRAM were performed 
using wobbler switching and a dual-polarization mode. At 
APEX position-switching was used with the reference position 
located at +2' in azimuth. The JCMT observations were also per- 
formed using position-switching with a 2' throw. Regular point- 
ing checks were performed during all observations and typically 
found to be consistent within a!3"of the pointing model. 

The data was reduced using CLASS, Stariink, and XSQ 
by subtracting a first order polynomial baseline fitted to the 
emission-free channels, and then binned to improve the signal- 
to-noise ratio. The typical spectral resolution of the reduced 
data is Ikms"'. The raw spectra are stored in the T* scale, 
where is the antenna temperature corrected for the atmo- 
spheric attenuation using the chopper-wheel method. The inten- 
sity scale is subsequently given in main-beam brightness temper- 
ature, - Ty-rfmb, where rf^b is the main-beam efficiency. The 
adopted main-beam efficiencies are given in Table |2] together 
with frequency, the energy of the upper level, and main beam 
FWHM for the respective transition. The uncertainty in the ab- 
solute intensity scale is estimated to be about +20%. Velocities 
are given with respect to the Local Standard of Rest (LSR). 



literature. When examining the previously published data, we 
found a large scatter in the reported CO(7 = 2— > 1) intensities 
and consequently these data were not used in the final analy- 
sis. The new CO( 7 - 2^ 1) data is gene rally stronger than the 
data discarded in iRamstedt et al.l (|2006|) . As an example, the 
new CO(7 = 2— > 1) line intensities obtained at IRAM are, on 
average, a factor of about three stronger than those reported in 
Sahai & Liechti (1995 ). 

The observed line profiles are in general of good quality 
(S/N ~ 5 or larger) and can, to a first approximation, be recon- 
ciled with those expected from a smooth spherically symmetric 
outflow. There is no evidence for any "detached-shell" sources in 
the sample (see e.g. the carbon star sample results of Olofsson 
et al. 1993). A closer inspection shows that some of the stars 
(e.g., W And, W Aql, DY Gem, and R Gem) have line profiles 
with weak wings extending beyond the parabolic line profile. 
This might be indicative of recent mass-loss-rate modulations in 
these sources or of asymmetric outflows. The six observed lines 
toward T Cet all show an asymmetry with stronger emission on 
the red-shifted side, similar to what is observed due to strong 
self-absorption ( iHuggins & HealvllI986l) . However, the star has 
a derived mass-loss rate of only 4x10"^ M© yr ' and all lines are 
found to be optically thin in the model. Consequently, the model 
does not reproduce the observed line profiles since strong self- 
absorption requires optically thick lines. Most likely, this is an 
indication that T Cet has an asymmetric outflow. 



3.3. SiO line profiles 

A total of 26 stars were detected in circumstellar SiO line emis- 
sion. All new spectra are presented in Appendix iBl and the inte- 
grated intensities and peak main-beam brightness temperatures 
are given in Tables |RT]|B3] The observed spectra are generally 
of good quality (S/N ~ 3 or lar ger). As previosly discussed by 
iGonzalez Delgado et al.l ( |2003|) . the SiO line profiles are often 
found to be narrower than the CO line profiles. However, many 
of the SiO lines show weak wings extending beyond the main 
emission feature and the total velocity widths are similar to those 
of the CO lines. It should be kept in mind that for poorer quality 
data, the uncertainty in the line width can be several km s ' . 



3.4. Dust continuum emission 

The SEDs were constructed using J, H, and K band data 
(2MASS) and IRAS fluxes. 



3.2. CO line profiles 

All new spectra are presented in Appendix [A] and the inte- 
grated intensities and peak main-beam brightness temperatures 
are given in Tables IaTIIA. 31 Ou r previous analysis of CO emis- 
sion lines from S-stars (Ram stedt et alj 12006) was based on 
new APEX observations in combination with data from the 



The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 

^ XS is a package developed by R Bergman to reduce and analyze 
a large number of single-dish spectra. It is publicly available from 
[ftp : //yggdrasil . oso . Chalmers . se , 



4. Radiative transfer modelling 

4. 1 . The circumstellar model 

The CSE is assumed to be spherically symmetric and formed 
by a constant mass-loss rate. It is assumed to be expanding at a 
constant velocity, derived from fitting the CO line widths, and 
to have a micro-turbulent velocity distribution with a Doppler 
width of 1 .0 km s ' . In addition, a thermal contribution to the lo- 
cal line width is added, based on the derived kinetic temperature 
of the gas. The density structure is obtained from the conserva- 
tion of mass. 
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4.2. Dust emission modelling 

The dust radiative transfer was solved using the publicly ay ail- 
able code DUST^Q. Amor phous carbon g rains (Suh 200C|) or 
amorphous silicate grains (iJusttanont & Tielen s 1992) were as- 
sumed based on the IRAS low-resolution spectra (LRS) classifi- 
cation according to lVoUc & Cohenl ( ll989l) . Three parameters are 
adjustable when fitting the model SEDs to the observed flux den- 
sities; the dust optical depth at 10 fim, tio, the dust temperature 
at the inner radius of the dust envelope, TdCn), and the stellar 
temperature, T*. One large grid per dust type was calculated. 
The optical depth at lOj-im ranges from 0.01 to 3.0 in steps of 
10 %, the inner dust temperature is varied from 500 to 1500 K, 
and the stellar temperature from 1800 to 2400 K, both in steps of 
lOOK. Once the grid is calculated, the solution can be scaled to 
the luminosity and distance of any star and a fit to the observed 
flux densities can be found ( llvezic & Elitzuill997h . For simplic- 
ity, the dust grains are assumed to be of the same size with a 
radius of 0. 1 yum, and have a density of 2 g cm""* (carbon grains) 
or 3 gem""* (silicate grains). 

4.3. Physical properties of the gas content in the CSEs 

The adopted method for the radiative transfer analysis of the cir- 
cumstellar CO radi o line emission has been de scribed in detail 
in p revious articles |Schoier & OlofssonI (1200 ih for carbon stars, 
and lOlofsson et al.](l2002l) forM-tvpe stars] and therefore only a 
short description is given here. 

The non-LTE radiative transfer code is based on the Monte 
Carlo method. The code has been benchmarked to high accu- 
racy a gainst a wide variety of molecular-line radiative transfer 
codes Ivan Zadelhoff et alj2002tlvan der Tak et al.ll2007h . In the 
excitation analysis of the CO molecule, 41 rotational levels are 
included in each of the two lowest vibrational states (v=0 and 
v=l). Energy levels and radiative transition probabilities, as well 
as collisional rate coefficients for collisions between CO and H2, 
are taken from^choier et al. (200513. When weighting together 
collisional rate coefficients for CO in collisions with ortho- and 
para-H2 an ortho-to-para ratio of 3 is adopted. Radiation from 
the central source (assumed to be a blackbody) and the cosmic 
microwave background is included. When the dust optical depth 
could be constrained in the dust radiative transfer analysis (see 
Sect. 15. II ). i.e., there is a significant amount of themal dust grains 
present in the wind, the corresponding modification to the radi- 
ation field is also included in the molecular excitation analysis. 

The energy balance equation for the gas is solved self- 
consistently. Cooling is due to line coohng from CO and H2 
and the adiabatic expansion of the gas. The dominant gas- 
heating mechanism is collisions between the H2 molecules 
and the dust grains. Photoelectric heating is also included and 
it is mostly important in the outer parts of the CSE. The 
free parameters describing the dust; the dust-to-gas mass-loss- 
rate ratio, *P, the average density of an individual dust grain, 
Pg, and the du st grain radius, a^, are co mbined in the pa- 
rameter h [see ISchoier & OlofssonI (1200 lb for a definition]. 
The /j-parameter thus gives a measure of how efficiently the 
gas is heated due t o coll isions with dust grains. Following 
ISchoier & OlofssonI ( 1200 lb . we have adopted an average effi- 
ciency factor for momentum transfer, which is constant through- 
out the CSE, <e,-p)"i =3 X 10-2. 



* http://www.pa.uky.edu/~moshe/dusty/ 
^ http://www.strw.leidenuniv.nl/~moldata 



The radial CO abundance distribution is estimated using the 
model presented in Mamon et al. (1988). The initial abundance 
of CO relative to H2, is assumed to be 6 x 10 The inner radius 
of the CSE was taken from the SED fitting when the model could 
be constrained, and is otherwise assumed to be 5 r*. A change in 
the inner radius with a factor of two will afi'ect the resu lting line 
intensities by less than 15% ( ISchoier & Olofssonll2001b . 

The mass-loss rate and the /z-parameter are the remaining 
free parameters in the CO Une modelling. 

4.4. SiO line radiative transfer modelling 

The same non-LTE Monte Carlo radiative transfer code is used 
for the analysis of the SiO radio line emission, where 41 ro- 
tational levels are included in each of the two lowest vibra- 
tional states (v=0 and v=l). Energy levels and radiative transi- 
tion probabilities, as well as collisional rat e coefficients for col - 
lisions between SiO and H2, are taken from lSchoier et alj(l2005b . 
The physical properties of the CSEs (mass-loss rates and radial 
gas temperature distributions) are taken from the results of the 
CO modelling. All other assumptions are the same. Also here, 
a decrease in the inner radius by a factor of two will result in 
a small change (<10%) of the model line intensities. The only 
free parameters when fitting the SiO line emission is the SiO 
abundance at the inner radius, and, when several lines were ob- 
served, the outer radius of the SiO emitting region (Sect. 14. 5t . 
These two parameters can be used to test chemical models of the 
inner CSE and photodissociation models for the outer envelope, 
respectively. 

4.5. The radial SiO abundance distribution 

The radial SiO abundance distribution, i.e., the ratio of 
the number densities of SiO molecules to H2 molecules, 
/ = n(SiO)/n(H2), is assumed to be described by a Gaussian 









- 




f) 


v 







For 18 out of 26 detected stars, too few transitions are observed 
in order to constrain both the abundance and the size of the emit- 
ting region. For these stars the extent of the SiO envelope is 
assumed to scale wi th the density of the CSE accor ding to the 
scaling law found by lGonzalez Delgado et al.l (l2003l) . 

(M\ 

log re = 19.2 + 0.48 log I — I , (2) 

where M is the mass-loss rate and Ve the gas expansion velocity 
of the wind, both found from the CO line modellin g. For a dis- 
cussio n about the val i dity of Eq. |2| see Gonzal ez Delgado et al.l 
(120031) . ISchoier et all (l2006bl) . and Sect. EH 

5. Results 

5.1. Dust optical depth and radial dust temperature 
distribution 

As known since previous studies of S-stars o n the AGB (e.g., 
ISahai & Liechtill995l:lJorissen & KnappI 19981) . the dust content 
of S-star CSEs is low. The dust radiative transfer analysis per- 
formed here show that the dust optical depth can be constrained 
for only 12 out of 40 sample stars. For the other stars the SED is 
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dominated by the emission from the central star and DUSTY has 
been used to determine only the stellar effective temperature (in 
addition to the distance/luminosity of the source). The best-fit 
model is found by minimizing 



Table 4. Mass-loss rates, /z-parameters, and terminal gas ex- 
pansion velocities, Vg, derived from fitting the radiative transfer 
model to the CO line observations. 



= 2 



mod,/! 



■ Fohs.A) 



(3) 



where is the flux density and cr^ the uncertainty in the mea- 
sured flux density at wavelength A. The summation is done over 
all independent observations. The reduced for the best-fit 
model is given by 



2 _ ^ min 
^red ~ 



(4) 



where p (the number of adjustable parameters) is 3 in our case. 

Table[3]shows the LRS classification, the assumed dust type, 
the stellar radius, r*, the inner radius of the CSE, r^, and the 
stellar temperature, . The inner radius is assumed to be 5 r* for 
all stars where the dust optical depth could not be constrained. 
Table [3] also gives the temperature at the inner radius of the dust 
shell, rd(ri), the dust optical depth at lOyum, rio, the dust-to-gas 
mass-loss-rate ratio, ^', and the dust mass-loss rate, Md, for the 
stars where this could be calculated. The dust mass-loss rate is 
given by 



In 

Md - Vd,oc n Tdy, 



(5) 



where x%r,v is the grain cross section per unit mass at a given 
frequency v, Vd,oo is the dust expansion velocity, and Tdv is the 
dust optical depth at a given frequency. A Xg^AQiim of 1.1x10"' 
and 3.5x10^ has been used for the carbon stars and M-type stars, 
respectively. The dust expansion velocity is calculated from the 
gas expansion velocity derived in the CO line modelling and the 
drift velocity between the dust and gas particles, Vdr, given by 



Vdr = 



Mc 

(lKwoklll975h . 



1/2 



(6) 



5.2. Mass-loss rates, radial gas temperature distribution and 
kinematics 

All results from the CO radiative transfer analysis, mass-loss 
rates, /i-parameters, gas expansion velocities, stellar velocities, 
vlsr, X^ed ^'^^ '-^^ best-fit models, and the number of observa- 
tional constraints, A^, are given in Table ID 

For slightly more than half of the stars, enough data (A^ > 3) 
is available to constrain both h and M. Models are calculated 
for a large number of mass-loss rates and /i-parameters and the 
best-fit model is found by minimizing 



X 



(^mod,i ^obs,i) 



(7) 



where /^od and /obs are the integrated intensities of the modelled 
and observed lines, respectively, cr is in most cases dominated 
by the calibration uncertainties assumed to be 20%. For obser- 
vations with poorer S/N, cr is set higher according to the quality 



Source 


M 

[lO-'Moyr-'] 


h 


[kms"'] 


vlsr 
[km s"'] 




N 


R And 


6.6 


0.4 


8.3 


-16.5 


1.3 


5 


W And 


1.7 


0.5 


6.0 


-35.0 


1.1 


5 


Z Ant 


0.9 


0.2 : 


6.0 


-16.0 


0.4 


2 


VX Aql 


1.0 


0.2 : 


7.0 


6.7 




1 


W Aql 


22 


1.0 


17.2 


-26.0 


0.6 


5 


AA Cam 


0.4 


0.2 : 


3.4 


-46.8 


2.8 


2 


T Cam 


1.0 


0.2 : 


3.8 


-11.7 


0.1 


2 


S Cas 


35 


0.2 


20.5 


-31.0 


3.6 


4 


V365 Cas 


0.3 


0.2 : 


6.2 


-2.1 




1 


WY Cas 


11 


0.5 


13.5 


7.0 


6.6 


3 


AM Cen 


1.5 


0.2 : 


3.2 


-27.0 


2.0 


2 


TTCen 


25 


0.5 : 


20.0 


4.0 


7.4 


2 


UYCen 


1.3 


0.2 : 


12.0 


-26.0 


0.2 


2 


V386 Cep 


2.0 


0.2 : 


16.0 


-51.8 




1 


TCet 


0.4 


0.2 


5.5 


22.0 


1.2 


6 


AACyg 


2.9 


0.2 


4.5 


27.5 


2.1 


3 


ADCyg 


2.1 


0.2 : 


8.0 


21.0 


0.1 


2 


RCyg 


6.3 


0.5 


9.0 


-18.5 


5.0 


5 


xCyg 


3.8 


0.2 


8.5 


9.0 


2.0 


6 


TVDra 


0.5 


0.2 


4.7 


22.0 


7.8 


3 


DYGem 


7.0 


0.2 


8.0 


-16.7 


2.0 


3 


RGem 


4.4 


0.8 


4.5 


-60.0 


3.7 


4 


ST Her 


1.3 


0.2 


8.5 


-4.5 


1.4 


4 


RXLac 


0.8 


0.2 : 


6.5 


-15.4 




1 


GILup 


5.5 


0.5 : 


10.0 


6.0 




1 


RLyn 


3.3 


0.2 


7.5 


16.0 


0.5 


5 


YLyn 


2.3 


1.0 


7.5 


-0.5 


3.9 


3 


SLyr 


20 


0.5 : 


13.0 


49.0 


7.6 


2 


FUMon 


2.7 


0.2 


2.8 


-42.0 


0.6 


3 


RZPeg 


4.6 


0.2 : 


12.6 


-23.4 




1 


RTSco 


4.5 


0.5 : 


11.0 


-47.0 


3.2 


2 


ST Sco 


1.5 


0.2 


5.5 


-4.5 


2.0 


3 


RZSgr 


30 


0.6 


9.0 


-31.0 


1.4 


3 


STSgr 


2.0 


0.5 : 


6.0 


55.7 


0.1 


2 


TSgr 


1.4 


0.2 


7.5 


9.0 


3.8 


3 


EP Vul 


2.0 


0.2 


5.7 


0.0 


8.9 


3 


DK Vul 


2.5 


0.6 


4.5 


-15.0 


5.7 


4 


AFGL 2425 


3.0 


0.2 


8.7 


57.0 


0.2 


3 


CSS2 41 


5.8 


0.2 


17.0 


21.5 


2.4 


3 


IRC-10401 


3.5 


0.2 : 


17.0 


19.0 


13.4 


2 



of the data (although never above 30%). The reduced;\^- is given 
by Eq.m where p -2. 

For the remaining stars, we have assumed the value of h 
depending on the stellar luminosity; h-Q.2 forLi_<_5000Lj7j_and 
h-0.5 for L j, > 5000 Lp [in accordance with lSchoier & OlofssonI 
(120011) and lOlofsson etal.l ( l2002l) l. This is indicated by a colon 
in TablelHand the reduced is given by Eq.|4]with p - I. 

Generally, the mass-loss rates are well constrained, and the 
observed CO lines are well-reproduced by the model. For many 
of the stars, it is not possible to find an upper limit to the h- 
parameter. This is not surprising, since the excitation in these 
low-mass-loss-rate stars is dominated by the radiation from the 
central star and is not very sensitive to the temperature structure 
of the gas and the adopted dust parameters. For some stars the 
X^^^ of the best-fit model is rather large. The derived is very 
sensitive to the line ratios, and any calibration errors in one line 
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(larger than our estimated uncertainties) might result in a large 
X^.^^ (see Fig. |2l CO(7 = 2^ 1)). Therefore it is difficult to de- 
termine whether inadequacies in the model or calibration errors 
are responsible for the poor fit in some cases. We find no system- 
atic trends with particular lines or telescopes. Figure|2]shows the 
best-fit model from the radiative transfer analysis overlayed on 
the observed line profiles for;^^ Cyg. The;^'--map is also shown 
with the innermost contour representing the Icr- level. This fit is 
good, but not perfect, e.g., in the details of the line shapes, and 

the main reason for this is not easily identified. 

Co mpared to the mass-loss rates found in iRamstedt et alj 
(I2OO6I) . the larger quantity of data makes the new mass-loss- 
rate estimates more reliable. In particular, it is possible to de- 
termine the /i-parameter for a larger number of stars. A com- 
parison shows that the inclusion of the new data in the analysis 
has changed the estimated mass-loss rates by less than a factor 
of two. For seven stars the change is slightly larger There is no 
systematic change upwards or downwards. 




5.3. SiO abundances and radial distribution 

All results from the SiO radiative transfer analysis, SiO frac- 
tional abundances, /o, the size of the SiO emitting regions (as 
derived from Eq. |2]i, r^, SiO gas expansion velocities, Ve(SiO), 
'^red ^'^^ '■^^ best-fit models, and the number of observational 
constraints, A^, are given in Table |5] Again, the best-fit model 
is found by minimizing Eq.|2]and;^'j:^j is defined by Eq.Hl where 
p=l when Te is given by Eq. [S] The derived abundances range 
from 4 X 10"^ to 1.4 x 10""*. The median value is 6 x 10"''. 
Figure[3]shows a histogram of the circumstellar SiO abundances 
for the S-type stars corripared to those found for M-type stars 
( Gonzalez Delgado et al.ll2003h and carbon stars dSchoier et alJ 



l2006bh . The spread of the SiO abundance distribution for the S 
type stars [as measured by the ratio between the 90th percentile 
and the 10th percentile] is about a factor of 25. 

The SiO gas expansion velocity is found to be, on average, 
20% sma ller than the CO gas expansion velocity (for 19 out of 
26 stars). [Gonzalez Delgado et alJ (l2003h discuss the possibility 



Fig. 3. Circumstellar SiO fractional abundances (/o) for the S- 
type (solid, g reen line; 26 stars), M-type (da shed-dotted, blue 
line; 45 stars; [Gonzalez Delgado et alj(2003|)). and carbon star 
(dashed, red line; 19 stars; dSchoier et al .I l2006bl) ) samples. 



that the SiO emission probes the very outer parts of the gas ac- 
celeration zone and that this would explain the discrepancy in 
the line width. SiO typically probes regions closer to the star by 
about a factor of 5-10 compared to CO. A trend of narrowing 
of the SiO lines with higher frequency might then also be ex- 
pected, since rotational transitions involving higher energy lev- 
els (Table |2]i tend to probe warmer and denser regions closer 
to the star. In addition, the line width can be reduced due to 
strong self-absorption on the blue-shifted side of the SiO lines 
that are optically thick. We do not find a trend with narrowing of 
higher frequency lines nor with the mass-loss rate [as would be 
expected if self-absorption was the predominant effect] and can 
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Fig. 4. The observed SiO line profiles toward R 
And and W Aql (histogram) overlayed with the 
results from the radiative transfer model (solid 
line) for the best-fit models with =0.6 and 
0.2, respectively. LSR velocity scales are used. 
The ;('^-maps of the full analysis for the two 
stars are also shown. 



therefore not firmly conclude that the SiO emission probes the 
acceleration zone. 



Table |6] shows the results (with Icr-errors) of the best-fit 
models when is left as a free parameter. In this case, the mod- 
elling generally gives a good fit to the observed SiO lines (see 
Fig.UJi withx^^^ on the order of unity. From the;^'^- analysis we 
find that the derived SiO abundances are generally determined 
within a factor of a; 4 (except for R And, Fig.|4]i. S Cas has a large 
;f^^^-value and the same can be observed for the result from the 
CO model (Tabled. This might indicate that S Cas is not well 
reproduced by a spherically symmetric wind model. The same 
can be noted for the SiO model ofx Cyg. However, both models 
are based on rather noisy SiO(y = 2 — »1) lines from OSO, and the 
line ratios between this Une and the other available lines are not 
reproduced. Thus, we can not decide whether the bad fit is due 
to the model being incorrect for these stars, or observational un- 
certainties. For T Cet, only the SiO(y = 8— »7) line was observed 
(at APEX). A fit to the width of the line gives Ve(SiO) = 16.0 km 
s Since only one line is available and the data is rather noisy 
(see Fig lB.lt , we have chosen to model the star with the expan- 
sion velocity derived from the CO model (N=6), Ve=5.5 km s"' . 
For lRC-10401 enough lines are available to try to fit the size of 
the SiO envelope and the abundance simultaneously. However, 



we are not able to constrain for this object and have chosen to 
use Te from Eq.|2l 

To be able to compare our res ults to what has previo usly been 
found for M-type stars (Gonzale z Delgado et al.ll20()3[) and car- 
bon stars dSchoier et al.l2006b) we have decided to derive all SiO 
abundances using Eq.|2]for the size of the SiO emitting region in 
order to be consistent. Figure|5]shows a comparison between the 
radii derived when is left as a free parameter and those using 
Eq. |2] The error bars correspond to the Icr limits. Only results 
from models with SiX^^d < 2-5 and where the SiO(7 = 8^7) line 
was available are shown in Fig.|5] We conclude that our results 
are consistent with Eq. |2] and that the size of the SiO emitting 
region does not differ significantly between the dif ferent chemi- 
cal type s. This is in agreement with the results of ISchoier et al.l 
(l2006bl) . 



6. Discussion 

6. 1 . Mass loss of S-type AGB stars 

Figure [6^ shows our mass-loss-rate results for the S-type 
AGB star sample compared with previous results for car- 
bon stars (Schoier & O lofsson 2001) an d M-type AGB stars 
jOlofsson et al.1 120021: iGonzalez Delgado eTaDllool). We con- 
firm here the result of [Ramstedt et al., (i2006h that the mass-loss- 



S. Ramstedt et al.: Circumstellar molecular line emission from S-type AGB stars: Mass-loss rates and SiO abundances 



9 



Table 5. Circumstellar SiO fractional abundance, /o, radius of 
the SiO line-emitting region from Eq. |2l r^, and gas expan- 
sion velocities derived from fitting the SiO line radiative transfer 
model to the observations, \'e(SiO). The reduced for the best- 
fit model and the number of observational constraints, A^, are 
also given. 



Source 
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/I 
4 


I V ura 


Zu 


Z.4 


/I ^ 
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O.J 


O.J 
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1.4 
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4.0 
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3 
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2 


STSco 


1.0 


3.7 


5.5 


7.9 
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17.0 
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rate distributions for the three different samples are very simi- 
lar. The median value for the S-type stars using the new data is 
2.7x10"^ Mq yr"', and the M-type and carbon stars both have a 
median of 3.0x10"^ Moyr"'. There are possibly fewer S-type 
stars with high mass-loss rates and due to our sample selec- 
tion criteria (see Sect.|2]i S-type stars with no or very little mass 
loss will be missed. The derived mass-loss rates depend on the 
adopted CO fractional abundance, which diff'ers for the three 
chemical types in accordance with stellar atmosphere models. 
Figure |6j5 shows the gas expansion velocity distribution derived 
from the CO line widths for the three samples. The median gas 
expansion velocity for the S-type stars is 8 kms"', and 7.5 and 
11 kms"^ for the M-type and carbon stars, respectively, indicat- 
ing that the carbon stars have CSEs with higher gas expansion 
velocities. Finally, Fig.|6}; shows the relation between the mass- 
loss rates and the expansion velocities for the three samples. We 
find no apparent difference depending on chemistry (apart from 
the higher gas expansion velocities in the carbon star sample) 
and suggest that this points to a mass loss that is driven by the 
same mechanism(s) in the S-type, M-type and carbon AGB stars. 

Another indication of this can be seen in Fig. [T] where the 
mass-loss rate (a and b) and expansion velocity (c and d) are 
plotted ag ainst the stellar var i ability period. As akeady dis- 
cussed by ISchoier & OlofssonI ( 1200 lb for the carbon star sam- 
ple, the mass-loss rate increases with the period of the star for 
all chemistries. The same trend can be observed for the expan- 
sion velocity. If the efficiency for driving a wind by radiation 



Table 6. Circumstellar SiO fractional abundance, /o, and the ra- 
dius of the SiO line-emitting region, rg, with Icr-errors, from 
the best-fit model. The reduced for the best-fit model and the 
number of observational constraints, A^, are also given. 
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Fig. 5. Circumstellar SiO envelope sizes (rg) estimated from 
the best-fit models plotted against the wind density (M/vg). 
The solid line is derived from Eq. |2l which is based on a fit 
to SiO envelope estimates for a sample of M-type AGB stars 
(iGonzalez Delgado et al.ll2003l) . The dashed line shows the pre- 
diction using a simple photochemical model (see Sect. 16.3b . 



pressure on dust was weaker in the S-type stars, due to, e.g., a 
low dust-to-gas ratio or a dust type incapable of driving a wind, 
the S-type stars would occupy a different parameter space than 
the M-type stars and the carbon stars in these plots. However, it 
is not possible to distinguish between the three different types. 

Q.2. Dust-to-gas ratio 

The average of the derived dust-to-gas mass-loss-rate ratios is 
about 2.8 X 10"^ for the S-type AGB stars in Table[3][V386 Cep 
excluded; this is likely the same type as GX Mon discussed 
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in iRamstedt etaLl (l2008l) 1. iGroenewegen et all (1 19991) derived 
dust-to-gas rati os of 48 M-type stars and found an average of 
5.8x10"^, and IRamstedt et all ( l2008h derived dust-to-gas ratios 
of four high-mass-loss-rate M-type stars resulting in, on average, 
2.8 X 10"^. For the carbon stars th e average dust-to-gas ratio s 
might be slightly lower [2.5 x 1 0~^, IGroenewegen et al] (Il998l) . 
2 x 10-\ IRamstedt et all (l2008l) 1. We conclude that the appar- 
ently low dust content in S-type AGB stars reflects the fact that 
they, for the most part, are low mass-loss rate objects. Their dust- 
to-gas ratios seem to be in agreement with what is derived for 
AGB stars of other chemical types, implying that the dust for- 
mation eflSciency is similar in all three chemical types. 



6.3. Circumstellar SiO abundances and constraints on 
chemical models 

The SiO abundances derived for the sample of S-type AGB stars 
range between 4x10"' and 1 .4 x 10""*. All stars (except RZ Sgr, 
which has the highest wind density) have abundances above or 
well above the expected therm al equilibrium value for C/O - 1 
(6.4 X 1 0"'. ICherchneg (l2006l) ). The median value of 6 x 10"* is 
almost an order o f magnitude larger than the equilibrium value. 
ICherchneffl(l2006h derives an SiO abundance at 5 of 3.3 x 10 
for C/O = 1 in a non-equilibrium shock-chemistry model. The 
SiO abundance derived in a shock-chemistry model is sensitive 
to the specific parameters of the star, such as the pulsational 
period and the shock velocity. The periods of the S-type AGB 



stars in our sample range between about 100-600 days but 
we find no obvious correlation with the derived SiO abundance. 
However, in a real, dynamic stellar atmosphere the chemistry is 
most likely much more complicated than in the present chemical 
models, and simple dependencies on single parameters, like the 
period, may not exist. We find no coiTelation between the LRS 
class and the derived SiO abundance for the S-type stars. 

The derived circumstellar SiO abundances for the sam- 
ple of S-type AGB stars are compared to those obtained 



for carbon stars ([S choier et al. 2006b) and M-type AGB stars 
(iGonzalez Delgado et al. 2003) in Fig. [3] The three distributions 
are very similar suggesting that the abundance of SiO in a cir- 
cumstellar chemistry is not very sensitive to the C/O ratio. In 
Fig.|8]the estimated abundances are plotted against a measure of 
the density of the wind (M/vg). For a specific density the SiO 
abundance can vary by up to two orders of magnitude. From the 
;if^-analysis we find that the abundances derived when also the 
size of the emitting region is a free parameter, are determined 
within a factor of ^ 4. Eq. |2] is derived using a relatively large 
number of stars with a large range in mass-loss rates and might 
give a statistically more accurate abundance. Similar investiga- 
tions of abundances of circumstellar molecules, for which the 
results are more in line with what would be expected in equi- 
librium chemistry (e.g., SiS and HCN), sho w a smaller range in 
the derived abundances fo r different stars dSchoier et al.ll2007l: 
ISchoier & Olofssonll2008l) . 
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For AGB stars with low density envelopes, such as the ma- 
jority of the stars discussed here, condensation of SiO molecules 
onto dust grains is not very effective and most probably not the 
cause for the observed scatter in the derived SiO abundances (see 
Sect. I6.5I 1. At chemical equilibrium a spread in the abundances 
is also to be expected depending on the C/O ratio and the tem- 
perature of the star. Around C/O - 1, the change in the amount 
of SiO formed is rather drastic, and between 0.95 and 1.05 the 
amount of SiO c an change with up to two orders of magnitude 
(lMarkwickl2000l) . If the temperature is varied between 2200 and 
2600 K for a given C/O-ratio, the SiO abundance will change 
less than an order of magnitude^. For the S-type stars alone it is 
not possible to conclude whether the spread in the derived abun- 
dances is indicative of a non-equilibrium chemistry or due to a 
spread in the C/O-ratio around 1 . However, given the results for 
all three chemical types, we conclude that the spread in the de- 
rived circumstellar SiO abundances most probably is real and 
indicative of a shock-chemistry in the formation zone of SiO. 

The derived SiO envelope sizes (see Table |6] and Fig.|5]l can 
be used to test photochemical models of the circumst ellar en- 
velop e . Using a simple photochemic al model CLindq vist et al.l 
I2OOOI: iGonzalez Delgado et al.l I2003L and references therein) 
with typical stellar, circumstellar and dust parameters for our 
sample of S-type AGB stars, and ado pting an unshielded pho- 
todiss ociation rate of 2.5xl0"^s ' ( iGonzalez Delgado et alj 
l2003h . we obtain the relation between SiO envelope size (re) 
and density measure (M/v'e) shown in Fig.|5](dashed line). These 
predictions agree well with the observed values and the relation 
described by Eq.|2] (solid line in Fig.|5]l used in the abundance 
estimates. 

6.4. SiO maser emission 

There are several indications that SiO maser emission orig inates 
from close to the stellar photosphere dReid & Moranl l 1981h . The 
energy of the masing state is about 1800K (v=l) or higher, 
and VLBA observations confirm that the masers are formed 
within a few stellar radii of the star (see for instance [Cotton et all 
I2OO6I and references there in). This corresponds to the region 
where the dust is for med jSalpeteij [T974t iDanchi et al.l Il994t 
iReid &Mentenlll997l) . Observations of SiO masers can there- 
fore provide information on whether there is SiO present in the 
dust formation region, and hence whether it can be one of the 
constituents of the dust formed. The non-detection of SiO maser 
emission toward carbon stars is thought to be indicating that the 
SiO molecules are formed further out in the wind than in M-type 
AGB stars. 

SiO maser emission has been searched for in 30 of 
our sample stars, and h as previously been detected in nine 
dJorissen & Knappll998l) . The results of our SiO(7 = 2^ 1, v=l) 
observations are given in Table IB. 31 and the detected lines are 
shown in Fig. IB. 31 In four stars the SiO maser has previously not 
been detected, and in another four stars we can not confirm pre- 
vious detections (R Cyg, R Lyn, Y Lyn and EP Vul), most likely 
a result of time variability. In Table IB73] we give distance inde- 
pendent upper limits for the sources observed. Unfortunately the 
reason for the non-detections is not possible to identify but time 
variability may play a role. 

Out of the 12 S-type stars in our sample classified as showing 
silicate emission in their LRS spectra, eight were observed and 
six were detected. Nine stars are classified as being featureless 
and out of these, six were observed and one was detected. Out 
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of the three stars classified as not showing any IR excess, two 
were observed and none was detected, and finally, out of the 13 
unclassified stars, seven were observed and two were detected 
in the SiO(/ = 2— > 1, v=l) maser line. These results are in line 
with a situation where the silicate stars have a C/O-ratio slightly 
less than one, and hence have enough SiO in their photospheres 
to produce SiO masers, while the others, in general, have C/O- 
ratios slightly larger than one, and hence too low photospheric 
SiO abundances to produce SiO masers. 



6.5. SiO cinemistry and implications for dust formation 

The formation of S iO in AGB winds with shocks is discussed by 
ICherchneS (l2006l) for different C/O-ratios. It is described how 
the formation of SiO is linked to the presence of OH and how 
the lack of OH in the inner wind of carbon stars explain their 
inability to form silicate dust despite their high SiO abundances 
further out in the CSE. 

For the S-type AGB stars, the situation is more complicated. 
The exact C/O-ratio is not known for the stars in our sample, 
but it is most likely a spread around 1 . The C/O-ratio of course 
has implications f or both the chemistry and the dust formed. 
ICherchne3 (l2006l) finds that the chemistry of the inner winds 
of stars with C/O - 0.98 is very similar to the M-type chemistry, 
while the chemistry of stars with C/O - 1.01 is very similar to 
that of carbon stars. In our comparison (Figs [3] and [8]) it is not 
possible to distinguish between the different chemical types in 
terms of their circumstellar SiO abundances. Given also the sim- 
ilarity of the mass-loss-rate distributions for the different chem- 
ical types, the most likely explanation is that the S-type stars in 
our sample have a spread in their C/O-ratio where some stars 
are more M-type like (forming silicate dust) and some are more 
carbon-star like (forming carbon dust). However, we are not able 
to exclude a scenario where two components of dust have formed 
in stars with a C/O-ratio close to unity. 

For M-type AGB stars ("Gonzalez Delgado et al1l2003h and 
carbon stars (Schoier et al. 2006b) there is a clear trend that 
the circumstellar SiO abundance gets lower as the density of 
the wind increases (Fig. [8]l, indicative of adsorption of SiO 
onto dust grains. For the sample of S-type AGB stars only an 
indication of such a trend may be seen, possibly due to the 
lack of high mass-loss rate objects. The dashed line in Fig. [8] 
shows a depletion curve based on a simple model presented in 
IGonzalez Delgado et al.l (l2003h . A relatively large scatter around 
the curve is expected since the condensation is sensitive to stel- 
lar, circumstellar, and dust characteristics. Further support for 
a deplet ion scenario comes fr om recent interferometric obser- 
vations. ISchoier et al.l (l2006ah modelled interferometric obser- 
vations of SiO(7 - 5^ 4) emission and infrared observations of 
rovibrational transitions of the extreme carbon star IRC+10216, 
and found that a two-component radial abundance distribution 
(a compact high abundance, pre-condensation, component com- 
bined with a more extended low abundance, post-condensation, 
component) is needed in order to explain the observations. 
Models of interfero metric observations of M-type AGB stars 
dSchoier et al.l l2004l) also show that a two-component radial 
abundance distribution gives a better fit to the observed data. 
The observations at hand for the S-type AGB stars do not allow 
a proper investigation of a scenario in which SiO molecules are 
incorporated into grains (see discussion in Schoieret al. (2007)). 
Instead, the abundances derived here should be considered as 
post-condensation abundances. 
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Fig. 8. Circumstellar SiO fractional abundances (/o) as a func- 
tion of wind dens ity (M/Ve) for S-type (gree n dots), M- 
type (blue squares) dGonzalez Delgado et al.ll20()3h . and carbon 
stars (red triangles) dSchoier et al.ll2006bh . The horizontal lines 
mark the abundan ces predicted from equilibrium chemistries 
jCherchneffl 12006*) . The dashed line shows the expected post- 
condensation abundance, /(oo), (scaled to 3.0x10"^, roughly the 
expected fractional abundance at 5 for low mass-loss rates 
when C/0=1) from a m odel including adsorption of SiO onto 
dust grains (see lGonzilez Delg ado et al...2003u for details). 



7. Conclusions 

We have modelled multi-transitional CO radio line data from 40 
S-type AGB stars. We have also modelled multi-transitional SiO 
data of 26 stars from the same sample. The results are compared 
to similar sample s of M-type AGB sta rs ( Olofsson et al.ll2002b 
and carbon stars dSchoier & Olofssonll2001h . We arrive at the 
following conclusions: 



- We find that the mass-loss rate distributions are very similar 
for the three chemical types (Fig.|6h), and so are the relations 
between mass-loss rate and expansion velocity of the stellar 
wind (Fig. [6]:). Further, it is not possible to distinguish be- 
tween the different chemical types when examining how the 
mass-loss rate correlates with the pulsation period of the star 
(Fig Q- The most likely explanation for the observed trends 
is a mass loss that is driven by the same mechanism(s) in 
all three chemical types. The only apparent difference is that 
the carbon star CSEs have higher gas expansion velocities 
(on average), most likely an effect of a higher acceleration 
efficiency of the carbon grains. 

- The derived dust-to-gas mass-loss-rate ratios for 1 1 of the 
sample stars show an agreement with what has previously 
been found for M-type AGB stars and carbon stars, imply- 
ing that the dust formation efficiency is similar in all three 
chemical types. 

- The median value of the estimated circumstellar SiO frac- 
tional abundances in S-type AGB stars is almost more than 
an order of magnitude higher than predicted by thermal 
equilibrium chemistry, suggesting that shock-induced non- 



equilibrium chemical processes are important in regulating 
the chemistry in the inner wind. 

The derived circumstellar SiO abundances for the S-type 
AGB stars range from 4 x 10"^ to 1.4 x 10""*. For a specific 
wind density (M/vg) there is a scatter of almost two orders 
of magnitudes, in accordance with what has previously been 
found for the other chemical types. In terms of their distri- 
bution of SiO circumstellar abundances, it is not possible to 
distinguish between the three different chemical types, and 
we propose that, although there is a large scatter, the circum- 
stellar SiO abundance is independent of the C/O-ratio for a 
given mass-loss rate and expansion velocity. The formation 
efficiency of SiO in a shock chemistry will be sensitive to 
other specific parameters of the star, like pulsational period 
and shock velocity. We believe that the scatter obtained in 
the SiO abundance estimates is real and indicative of a shock 
chemistry in the outer atmosphere. 

Previous analysis of M-type AGB stars 

Gonzalez Delgado et al.l l2003h and carbon stars 
Schoier et al. I I2006bh show a clear trend in that the 



circumstellar SiO abundance decreases as the density of the 
stellar wind increases. The trend is indicative of adsorption 
of SiO onto dust grains. The same trend can be suspected 
for the sample of S-type AGB stars, however, the number 
of high-mass-loss-rate S-type AGB stars is low and no firm 
conclusion can be drawn. 

Despite the limitations due to the selection criteria of the 
samples compared, we believe that our conclusions apply to 
AGB stars in general. 
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Table 3. LRS classification, adopted dust type (amorphous silicate, AMS, or amorphous carbon, AMC) , stellar radius, r*, inner 
radius of the CSE, n, stellar temperature, T*, temperature at the inner boundary, TdCn), dust optical depth at lO/zm, tiq, drift 
velocity, Vdr, dust mass-loss rate, Md, and dust-to-gas mass-loss-rate ratio, *F. 
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2.6 


1.3 


2400 












RZPeg 


C 


AMC 


2.6 


1.3 


2400 












RTSco 


E 


AMS 


3.3 


1.7 


2400 












ST Sco 


E 


AMS 


2.6 


1.3 


2400 












RZSgr 


F 


AMC 


2.6 


1.3 


2400 












ST Sgr 


E 


AMS 


3.1 


1.6 


2400 












TSgr 


F 


AMC 


3.1 


1.6 


2400 












EP Vul 


F 


AMC 


2.6 


1.3 


2400 












DK Vul 


F 


AMC 


2.6 


1.3 


2400 












AFGL 2425 




AMS 


4.6 


2.9 


1800 


600 


0.10 


8.4 


1.4 


4.8 


CSS2 41 




AMS 


4.6 


4.1 


1800 


500 


0.10 


8.5 


3.0 


5.2 


IRC-10401 




AMS 


4.6 


4.1 


1800 


500 


0.07 


11 


2.3 


6.6 
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Appendix A: New CO spectra 
Appendix B: New SiO spectra 



Source 


Transition 


-'mb 




Telescope 








[Kkms~'] 


[K] 


or Reference 


R And 


CO(l 


-0) 


5.7 


0.44 


BL94 




CO(l 


-0) 


34.9 


2.4 


I 




CO(2 


-1) 


135.5 


9.6 


1 




CO(2 


- 1) 


32.0 


2.5 


J" 




CO(3 


-2) 


43.0 


3.3 


J" 


W And 


CO(l 


-0) 


1.4 


0.12 


BL94 




CO(l 


-0) 


3.8 


0.4 


O 




CO(l 


-0) 


9.5 


0.9 


1 




CO(2 


-1) 


55.3 


4.7 


1 




CO(3 


-2) 


7.0 


0.6 


Y95 


Z Ant 


CO(2 


- 1) 


0.7 


0.07 


JK98 




CO(3 


-2) 


2.0 


0.6 


A 


VX Aql 


CO(l 


-0) 


0.7 


0.05 


SL95 


W Aql 


CO(l 


-0) 


28.4 


1.3 


N92 




CO(2 


- 1) 


54.8 


2.1 


K98 




CO(3 


-2) 


93.2 


3.5 


K98 




CO(3 


-2) 


136.2 


5.0 


A 




CO(4 


-3) 


199.5 


7.4 


J" 


AA Cam 


CO(2 


-1) 


1.8 


0.06 


SL95 




CO(3 


- 2) 


0.8 


0.17 


J 


TCam 


CO(l 


-0) 


2.0 


0.3 


I 




CO(2 


-1) 


8.1 


1.2 


I 


S Cas 


CO(l 


-0) 


13.8 


0.4 


O 




CO(l 


-0) 


38.2 


1.1 


1 




CO(2 


- 1) 


157.3 


4.8 


1 




CO(3 


-2) 


31.0 


1.1 


J" 


V365 Cas 


CU(Z 


- 1) 


0.7 


0.07 


SL95 


WYCas 


CO(l 


-0) 


3.5 


0.12 


O 




CO(l 


-0) 


12.8 


0.6 


1 




CO(2 


-1) 


64.4 


2.9 


I 


AM Cen 


CO(l 


-0) 


0.6 


0.08 


SL95 




CO(3 


-2) 


1.1 


0.24 


A 


TT Cen 


CO(l 


-0) 


1.9 


0.05 


SL95 




CO(3 


-2) 


13.9 


0.5 


A 


UY Cen 


CO(2 


- 1) 


1.0 


0.07 


SL95 




CO(3 


-2) 


1.6 


0.1 


A 


V386 Cep 


CO(2 


-1) 


8.4 


0.3 


I 


TCet 


CO(2 


-1) 


2.0 


0.28 


S 




CO(2 


-1) 


3.7 


0.4 


J" 




CO(3 


-2) 


6.1 


0.7 


A 




CO(3 


-2) 


7.2 


0.8 


J" 




CO(4 


-3) 


9.5 


0.97 


J" 




CO(6 


-5) 


10.0 


1.1 


J" 



Table A.l. Observational results of circumstellar CO radio line 
emission. 
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Table A.2. Observational results of circumstellar CO radio line 
emission (continued). 



Table A.3. Observational results of circumstellar CO radio line 
emission (continued). 



Source 


Transition 


Imb 




Telescope 






[Kkms-'] 


[K] 


or reference 


AACyg 


C0(1 - 0) 


1.5 


0.11 


BL94 




C0(1 - 0) 


5.7 


0.7 


1 




C0(2 - 1) 


11 A 


2.3 


I 


ADCyg 


C0(1 - 0) 


0.9 


0.08 


1 




C0(2 - 1) 


3.8 


0.31 


1 


RCyg 


C0(1 -0) 


2.5 


0.15 


BL94 




C0(1 -0) 


4.4 


0.3 


O 




CO(l - 0) 


16.4 


1.0 


1 




CO(2- 1) 


70.2 


4.2 


1 




CO(3 - 2) 


14.6 


2.3 


S95 


y Cv2 

A. D 


CO(l - 0) 


27.2 


1.8 







CO(l - 0) 


70.9 


6.0 


1 




CO(2- 1) 


60.2 


4.2 


J" 




CO(2 - 1) 


265.8 


19.8 


1 




C0(3 - 2) 


119.5 


8.3 


J 




C0(4 - 3) 


134.6 


9.6 


J" 


TV Dra 


C0(1 -0) 


0.6 


0.02 







C0(1 - 0) 


1.5 


0.16 


1 




C0(2 - 1) 


4.1 


0.6 


I 


DYGem 


C0(1 - 0) 


4.4 


0.3 


1 




C0(2 - 1) 


17.1 


1.3 


1 




CO(3 - 2) 


2.6 


0.35 


A 


RGem 


C0(1 - 0) 


2.4 


0.3 







CO(l - 0) 


9.7 


1.3 


I 




CO(2- 1) 


34.1 


5.2 


1 




C0(3 - 2) 


6.5 


0.97 


A 


ST Her 


C0(1 - 0) 


1.5 


0.2 


O 




CO(l - 0) 


4.9 


0.24 


KJ94 




CO(2- 1) 


17.4 


1.05 


KJ94 




CO(3 - 2) 


5.1 


0.34 


K98 


RXLac 


C0(1 -0) 


2.0 


0.3 


GJ98 


GILup 


CO(3 - 2) 


4.6 


0.3 


A 


RLyn 


C0(1 - 0) 


0.3 


0.03 


BL94 




CO(l - 0) 


0.8 


0.1 







CO(l - 0) 


2.0 


0.18 


1 




C0(2 - 1) 


5.8 


0.5 


1 






n Q 


U.Uo 




YLyn 


C0(1 -0) 


4.1 


0.3 


o 




C0(1 -0) 


14.3 


1.1 


1 




C0(2 - 1) 


56.4 


4.4 


1 


SLyr 


C0(1 -0) 


1.9 


0.1 







C0(3 - 2) 


6.9 


0.32 


A 


FUMon 


C0(1 - 0) 


3.3 


0.7 


1 




CO(2- 1) 


7.0 


2.0 


1 




C0(3 - 2) 


1.3 


0.3 


A 


RZPeg 


C0(2 - 1) 


3.5 


0.21 


SL95 



Source 


Transition 




Tmb 


Telescope 






[Kkms-'] 


[K] 


or reference 


RTSco 


C0(1 -0) 


3.8 


0.20 


SL95 




C0(3 - 2) 


32.6 


1.8 


A 


ST Sco 


CO(l -0) 


0.6 


0.06 


BL94 




C0(1 - 0) 


1.4 


0.14 


SL95 




CO(3 - 2) 


4.5 


0.5 


A 


RZSgr 


C0(1 -0) 


6.4 


0.38 


S 




C0(2 - 1) 


10.7 


0.87 


S 




C0(3 - 2) 


11.7 


1.0 


A 


STSgr 


C0(1 -0) 


0.8 


0.66 


SL95 




C0(3 - 2) 


4.9 


1.0 


A 


TSgr 


C0(1 - 0) 


1.5 


0.17 


1 




C0(2 - 1) 


5.2 


0.4 


1 




C0(3 - 2) 


2.5 


0.2 


A 


EP Vul 


C0(1 - 0) 


0.9 


0.12 


BL94 




C0(1 -0) 


0.6 


0.08 


SL95 




C0(3 - 2) 


3.1 


0.35 


A 


DK Vul 


C0(1 - 0) 


1.3 


0.04 


BL94 




C0(1 -0) 


2.9 


0.6 


1 




C0(2 - 1) 


9.6 


1.7 


1 




C0(3 - 2) 


3.6 


0.7 


A 


AFGL 2425 


C0(1 -0) 


3.3 


0.2 


1 




C0(2 - 1) 


10.5 


0.75 


1 




C0(3 - 2) 


2.8 


0.25 


A 


CSS2 41 


C0(1 -0) 


2.6 : 


0.15 : 


1 




C0(2 - 1) 


4.7 : 


0.3 : 


1 




C0(3 - 2) 


2.5 : 


0.1 : 


J 


lRC-10401 


CO(l -0) 


2.9 


0.12 


GJ98 




C0(3 - 2) 


9.5 


0.45 


A 



A; APE X, 1; IRAM, J; JCMT, J°; J CMT arch ival data , O; OS O 
BL94: iBieging & Latteil ( fl994l) . Y95: lYound i fT99^. JK98; 
Jorissen & Knapd \l99m . SL95: ISahai & Liechtil ( 1 19951 ). N92; 



j 19981) . 

Nvmanetal. 1 |l992h. K9 8 



KJ94; 



Knapp etalJ ||1998^ 
Kahane & Jura Cl994i) , GJ98 



S95; Staneketal 
Groenewegen & de Jong 
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Table B.l. Observational results of circumstellar SiO radio line 
emission. 



Table B.2. Observational results of circumstellar SiO radio line 
emission (continued). 



Source 


Transition 


[K km s"' 


[K] 


Telescope 
or Reference 


Source 


Transition 


-'rab 

[K km s"' 


Tmb 

[K] 


Telescope 
or Reference 


R And 


SiO(2- 1) 


0.9 


0.06 


O 


RT Sco 


SiO(3 - 2) 


0.5 


0.04 


S/B98 




SiO(2- 1) 


2.5 


0.22 


1 




SiO(8 - 7) 


1.5 


0.19 


A 




SiO(5 - 4) 


13.4 


1.12 


1 














SiO(6 - 5) 


3.2 


0.31 


J 


ST Sco 


SiO(3 - 2) 


0.1 


0.02 


S/B98 




SiO(8 - 7) 


3.5 


0.39 


A 




SiO(8 - 7) 


0.5 


0.07 


A 


W And 


SiO(2- 1) 


1.7 


0.25 


1 


RZSgr 


SiO(3 - 2) 


0.3 


0.03 


S/B98 




SiO(5 - 4) 


6.8 


1.1 


I 




SiO(8 - 7) 


0.5 


0.10 


A 


W Aql 


SiO(2- 1) 


2.6 


0.08 


O 


STSgr 


SiO(8 - 7) 


0.6 


0.08 


A 




SiO(3 - 2) 


3.9 


0.15 


S/B98 














SiO(8 - 7) 


10.1 


0.52 


A 


TSgr 


SiO(8 - 7) 


0.5 


0.09 


A 




SiO(8 - 7) 


16.9 


0.84 


J 






















EP Vul 


SiO(2 - 1) 


0.4 


0.07 


I 


S Cas 


SiO(2- 1) 


1.7 


0.06 


O 




SiO(5 - 4) 


1.8 


0.27 


I 




SiO(2- 1) 


6.1 


0.22 


I 




SiO(8 - 7) 


0.7 


0.12 


A 




SiO(5 - 4) 


35.8 


1.5 


^ 






















DK Vul 


SiO(2 - 1) 


0.2 


0.03 




WYCas 


SiO(2- 1) 


0.6 


0.04 


1 




SiO(5 - 4) 


1.0 


0.14 


I 




SiO(5 - 4) 


3.3 


0.2 


I 














SiO(8 - 7) 


0.9 


0.11 


J 


AFGL 


SiO(2 - 1) 


0.9 


0.07 


I 












2425 










VJso Cep 


SiO(2 - 1) 


1.2 


0.07 






SiO(5 - 4) 


2.0 


0.20 






SiO(5 - 4) 


4.3 


0.26 






SiO(8 - 7) 


0.5 


0.08 


A 




SiO(8 - 7) 


2.1 


0.13 
























CSS2 41 


SiO(2 - 1) 


1.0 


0.05 




TCet 


SiO(8 - 7) 


0.7 


0.03 






SiO(5 - 4) 


2.4 


0.16 




AACyg 


SiO(2- 1) 


0.3 


0.05 




IRC- 


SiO(2 - 1) 


2.9 


0.14 






SiO(5 - 4) 


2.0 


0.3 




10401 






















SiO(5 - 4) 


2.2 


0.20 




RCyg 


SiO(2- 1) 


0.3 


0.03 






SiO(8 - 7) 


1.4 


0.12 


A 



A'Cyg 

TVDra 

DYGem 
RGem 

RXLac 

YLyn 

SLyr 



SiO(2- 1) 
SiO(5 - 4) 
SiO(8 - 7) 

SiO(2- 1) 
SiO(2- 1) 
SiO(5 - 4) 
SiO(8 - 7) 

SiO(2- 1) 
SiO(5 - 4) 

SiO(5 - 4) 

SiO(2- 1) 
SiO(3 - 2) 
SiO(5 - 4) 

SiO(2- 1) 
SiO(5 - 4) 
SiO(8 - 7) 

SiO(2- 1) 
SiO(5 - 4) 
SiO(8 - 7) 

SiO(2- 1) 
SiO(5 - 4) 



1.3 
6.3 
2.0 

4.5 
23.1 
72.3 
38.1 

0.6 
1.3 

1.1 

0.2 
0.1 
1.7 

0.3 
1.7 
0.6 

0.4 
2.6 
1.8 

0.5 
3.0 



0.14 
0.65 
0.20 

0.47 
1.94 

6.7 
3.56 

0.08 
0.21 

0.09 

0.06 
0.03 
0.30 

0.05 
0.30 
0.07 

0.07 
0.33 
0.14 

0.04 
0.19 



A; A PEX, 1; IRAM J; JC MT, O; OSO, S; SEST 
B98: lBiegingetal.l ( ll998h 



S/B98 
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Table B.3. Observational results of circumstellar SiO(y = 2-1, 
V = 1) emission. 



Source T^b cr r,„b X prev. 

[K] [K] [Kpc^]" det.'' 



R And 


0.52 


0.016 


4.7 X 


10* 


Y 


W And 


0.72 


0.027 


5.7 X 


10* 


Y 


VX Aql 




047 


O.O A 


lO"* 

iU 




W Aql 


0.17 


0.022 


9.0 X 


103 


Y 


AA Cam 




0.022 


4.2 X 


10* 




S Cas 


1 Q 

i.y 


075 


3.7 X 


10^ 


Y 






0.026 


3.0 X 


10* 




WY Cas 


0.15 


0.014 


5.4 X 


10* 






0.70 


0.019 


1.5 X 


IqS 




A A Tver 




07"^ 


1.6 X 


10* 


N 


An Cvo 




05^ 


1.5 X 


10^ 




K Lyg 




U.UUo 


H.J A 


lU 


V 


xCyg 


8.1 


0.082 


9.8 X 


10* 


Y 


TV Dra 


0.08 


0.012 


1.2 X 


10* 


N 


DYGem 




0.030 


4.2 X 


10* 




R Gem 


0.19 


0.015 


9.3 X 


10* 


N 


ST Her 




0.022 


6.0 X 


10' 




RXLac 




0.032 


9.3 X 


10' 




RLyn 




0.027 


6.0 X 


10* 


Y 


YLyn 




0.012 


2.4 X 


10' 


Y 


SLyr 




0.037 


1.6 X 


10= 




FUMon 




0.038 


6.9 X 


10* 


N 


EP Vul 




0.050 


3.9 X 


10* 


Y 


DK Vul 




0.043 


7.2 X 


10* 





" For t he non-detections we use 30" x D^. 
From lJorissen & Knapd ( 119981) and references therein. 
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Fig.A.l. New observations of circumstellar CO line emission. The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 
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Fig.A.2. New observations of circumstellar CO line emission. The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 
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Fig.A.3. New observations of circumstellar CO line emission. The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 
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Fig. B.l. New observations of circumstellar SiO line emission (v = 0). The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 



S. Ramstedt et al.: Circumstellar molecular line emission from S-type AGB stars: Mass-loss rates and SiO abundances 



23 




Fig. B.2. New observations of circumstellar SiO line emission (v = 0). The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 




Fig. B.3. New observations of circumstellar SiO line emission (v = 1). The source name is shown to the upper left, and the transition 
observed and the telescope used is shown to the upper right of each frame. LSR velocity scales are used. 



